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We develop a systematic classical framework to accommodate canonical quantization of geomet- 
ric and matter perturbations on a quantum homogeneous isotropic flat spacetime. The existing 
approach of standard cosmological perturbations is indeed proved to be good only up to first order 
in the inhomogeneities, and only if the background is treated classically. To consistently quantize 
the perturbations and the background, a new set of classical phase space variables is required. We 
show that, in a natural gauge, a set of such Dirac observables exists, and their algebra is of the 
canonical form. Finally, we compute the physical Hamiltonian that generates the dynamics of such 
observables with respect to the homogeneous part of a K-G "clock" field T. The results of this work 
£C) [ provide a good starting point to understanding and calculating effects that quantum cosmological 

spacetime in the background has on the quantum perturbations of the metric tensor and of matter 
fields. 
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I. INTRODUCTION 
A. Motivation: towards a new quantum theory 



The framework we present in this paper is classical, but the aim is introduction of quantum test fields on a quantum 
I • spacetime. While the theory of test quantum fields propagating on classical FRW spacetime is very well known, a 
£h ' theory of quantum fields on a quantum universe in the background is a scenario considered much less often. The first 
i , step in that direction was made in In that work, the quantum background was provided by the loop quantum 
cosmology model 1 of the homogenues isotropic universe characterized by the scale factor e a coupled to a homogeneus 
massless K-G field (playing the role of physical time), and the test field was a second K-G field 5(f>. This model 
*■ ■ was derived from the (suitably simplified) theory of two K-G fields, T(x) and <fi(x), coupled to gravitational field 
9fj,u(x): the authors started with the scalar constraint of the full theory, expanded it around the homogeneus solutions 
(with = 0), and dropped all the degrees of freedom except for (i) the scale a of the universe, (ii) the homogeneus 
j part of the first K-G field, and (iii) the perturbations S(f> of the second K-G field. All these three remaining 
degrees of freedom were coupled to each other and subject to a quantum scalar constraint defined by the truncated 
quantum constraint operator C — C a + C T (o) + Cs<t> — 0. (The vector constraint was satisfied automatically at the 
classical level in the test field approximation, therefore it was ignored at the quantum level.) That idea was later 
generalized to the Bianchi I quantum spacetimes @. 

The goal of that new approach to QFT on quantum spacetimes was to gain some insight into possible effects of the 
quantum nature of geometry on the propagation of test fields. A first intriguing conclusion came from the study of 
a mechanism of the emergence of a classical spacetime from such quantum system. The classical spacetime emerges 
as a metric tensor effectively felt by the modes of the quantum test field. It is obtained from the quantum dynamics 
of the modes, and differs from solutions to the classical Enstein equation with suitable quantum corrections. For 
that reason, it was later called "dressed" 10]. If the test K-G field is massless, then all its modes experience a single 
dressed metric, independently of their momenta. That phenomenon may be interpreted as the absence of so-called 
Lorenz symmetry violation. If the test K-G field is massive, on the other hand, then the dressed metric felt by each 
mode depends of the direction of its momentum. In particular, the dressed metric is not any longer space-isotropic, 
even in the case of an isotropic quantum universe in the background [111 ]. One might say that, from the point of 
view of each mode, the isotropy is broken by quantum geometry effects! These results were insensitive on possible 
choices of the quantum model for the homogeneous degrees of freedom, be them loop quantum cosmology |6|-l8j or 
Wheeler-deWitt quatum cosmology [l2l [l~3j . However, the findings outlined above were derived using a quite crude 
approximation, therefore one expects more from a systematic approach. 
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Since the concept of the test quantum field on the quantum cosmological background spacetime has proved to 
be quite fruitful, it is worth to extend it to perturbations of the gravitational field and to perturbations of the the 
scalar field Tpresent in the background spacetime. An attempt to achieve that goal was made in [13]. Therein, the 
approach of [l| is reconsidered. However, it is developed in a somewhat different direction. Namely, as a starting 
point a very well known standard perturbation theory of classical mathematical cosmology p^4l6| is taken. Each 
dynamical variable 7 (a coordinate in the phase space) is expanded as 

7 = 7 (0) +e57 (l) + i e 2^(2) + ... (1) 

where 7C) encodes the background part. The background part is decoupled from the higher order perturbations. 
Thus, the 1st order perturbations S^ 1 ' gain their own phase space and are subject to a theory introduced on a 
fixed background 7W. Such framework is powerful in cosmology and useful for example in the context of QFT on 
curved classical spacetime, because the background spacetime has been fixed from the beginning as a solution to the 
unperturbed Einstein equations, and the phase space of the system consists of the perturbation sector only. In [To| that 
framework is somehow quantized. In the consequence, one quantizes the perturbations order by order, decoupled from 
the background and from each other. What we want to consider instead, is a joint quantization of the perturbations 
and of the background (as is the case of [H, Q). Therefore, we will propose in the current paper to go in an alternative 
direction to that that of [l(|. To this end we need to develop a classical framework that keeps the original coupling 
between the perturbations and the background. Another proposal for a systematic development of the theory of 
quantum perturbations coupled to quantum background was made in |17| . The starting point of the authors is the 
classical framework of [HI available for the k = 1 cosmology (a spherical universe). Their framework combines the 
LQC quantization of the background with a unique quantization of scalar cosmological perturbations on a classical 
background with k = 1 [19j . That proposal is satisfactory from the point of view of deriving a theory of perturbations 
from the full theory. What we are looking for in the current paper is a similar classical starting point available for a 
flat universe. 



B. Our goal: a compatible classical framework 

The current paper 2 provides such classical framework, suitable for the project of studying QFT on quantum space- 
time. For the sake of self-sufficiency, we start from scratch. We address the full theory of gravity coupled to two K-G 
fields, and systematically develop a classical framework in which the perturbations and the background together set a 
phase space. To this end, throghout this paper we consider the full phase space of the system of the two Klein-Gordon 
fields coupled to the gravitational field. Thereon, we introduce a coordinate system which consists of coordinates 7' ' 
and coordinates £7. This coordinate system is defined globally on the phase space, and it is exact. At this point the 
reader should not be mislead by the symbol 5 in front of 7: nothing is small yet. In such phase space we consider 
the constraints of the full theory, reduce to the constraint surface, and there we study the gauge transformations 
generated by the constraints, which finally allow to restrict to the reduced physical phase space. We construct the 
Dirac observales, along with a 1-dimensional group of automorphisms of their algebra parametrized by the variable 
T'°), and find the generator of such group, i.e. the physical Hamiltonian. The structure that we needed for the 
quantization (and that we derive below) is the Poisson algebra of the Dirac observables and the physical Hamiltonian. 
Technically, our goal in this work is to find some basic Dirac observables which define the canonical Poisson algebra, 
and to express by them the physical Hamiltonian. 

Reading the outline of the programme we just gave, one could say that we made the dream of a canonical relavist 
come true! The catch is, in order to see explicit expressions we implement the canonical programme expanding 
everything in S^y and keeping track of the results of the calculations up to the 1st order only. Nevertheless, the 
framework allows to continue to the second order, and actually to any higher order. 



C. Comparison and contrast with the standard approach 

For the reader familiar with the standard approach to cosmological perturbations [l4l - fl6j , we outline here similarities 
and differences between the standard framework and our approach. 



2 The first ideas were presented in |20| . 
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As in the standard perturbative framework, we split every field into its background and perturbation part, 

7 = 7<°> + <y 7 (2) 

However, there is no higher order 5^ 2 \ since both 7 (°) and 5 7 are defined exactly on the full phase space of the system, 
and our perturbation part of the metric is not just the first approximation. More precisely: the expression ([2]) is not 
the first two terms of the Taylor expansion (there is no "small parameter" e, and both terms are finite), but rather an 
exact, unique decomposition giving rise to a certain coordinate system on the phase space. To give an example, for the 
K-G field T(x) sector, we will define to be the homogeneous part and ST(x) the inhomogeneous part: that split is 
always well defined, and does not involve any knowledge of the dynamics (contrary to Jl}, where given an exact T(x), 
the definition of is the part of T(x) that satisfies unperturbed K-G equation). Because of the kinematical nature 
of this decomposition, it follows that in our sense, 5T{x) does not involve any correction to the homogeneous part of 
the field: it is all absorbed in T^ ' (which in the consequence does not satisfy the unperturbed order K-G equation). 
Actually, if one considers the dynamics up to the first order, the variables 7^°- ) do satisfy the Klein-Gordon-Einstein 
equations (with <5 7 = 0), and the variables £7 satisfy the linearized equations on the background 7W. Thus, up to 
the 1st order, the two approaches agree. However, it will be clear that, if quadratic order were considered, then this 
would not be the case, since 7^ would know about "back-reaction of perturbations on the background" . 

The essential feature of our framework is that our y(°> and #7 are subject to the unchanged Poisson algebra of 
the full theory, contrary to j(°\...,j( n \... in standard cosmological perturbation theory. Starting with the correct 
classical Poisson algebra is relevant for the (future) canonical quantization. Of course, the price to pay is that now 
7^°) is not a fixed background: it is a dynamical background, its dynamics being generated by the same Hamiltonian 
that generates the dynamics for the perturbations <$ 7 . 

Another important difference between the two approaches is in the structure of the constraints and the gauge 
transformations. In our approach, the constraints and the gauge transformations are those of the canonical Klcin- 
Gordon-Einstein theory. On the other hand, in the standard perturbation theory, for the first order perturbations 
(57W only the linear part of the full constraints is considered a generator of gauge transformations. For that reason, 
we will see that the Mukanov-Sasaki variables Q and P - widely used in standard cosmological perturbation theory 
as the gauge-invariant observables - are no longer Dirac observables in our approach, even up to the 1st order. That 
is why we need to introduce our Dirac observables, show that their Poisson algebra is canonical, and endow it with 
dynamics derived explicitely up to the 1st order only. 

All these differences are unavoidable in our programme. In fact, this feature of treating the homogeneous and 
inhomogeneous parts on the same footing is precisely what motivates us into thinking that our proposed approach is 
more suitable for canonical quantization of perturbations and background. 



II. THE KLEIN-GORDON-EINSTEIN THEORY 

We consider the gravitational field coupled to two Klein-Gordon fields. The system is described by the following 
action: 



S = Id x\/—g 



- i/'S^ - V T (T) - \ g ^d^d v 4> - W) 



(3) 



2k 2* p " ^ ' 2' 
where k = 8ttG and R is the Ricci scalar of the gravitational field <? M „. We can distinguish three sectors: 

• the geometric (G) sector, associated to the metric g M „, 

• the time (T) sector, a Klein-Gordon "clock field" T, 

• the matter (M) sector, a Klein-Gordon field 4>. 

The field T is referred to as time, because the value of its spatially homogenous part (with respect to fixed coordinates), 
T™, will be used to parametrize a one dimensional group of automorphisms acting on the Dirac observables. 

To proceed with the canonical quantization of the theory, it is convenient to pass to the Hamiltonian formalism. 
The usual way to do so is to write the metric in the ADM form: 

g^dx^dx" — —(A 2 — q ab N a N b )dt 2 + 2q ab N b dtdx a + q ab dx a dx b (4) 

Here, q ab is the spatial metric, i.e. the metric that g M „ induces on the spatial (Cauchy) surface X C M . N and 
N a are called respectively the lapse function and the shift vector field, and characterize the spacetime geometry in 
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the directions transversal to £ embedded in M. In our paper, the latin indices run through the set {1, 2, 3} and are 
characteristic to the objects living in the tensor bundle of S. They are raised and lowered by using the spatial metric. 
In matrix form, the metric and its inverse are given by 



I -N 2 + N a N a N a 
N a q ab 

V 



-1/N 2 
N a /N 2 



„ab 



N a /N 2 
- N a N b /N 2 



(5) 



From here, we can perform the canonical analysis of the action ((3]). First, we plug ((4]) in ([3]) , then use the Codazzi 
equation to write R in terms of the 3-dimensional Ricci tensor and of the extrinsic curvature 



K ab = --{C n q)ab = -qjt [<i°-b - (£-N<l)ab\ = 
= ^ [-qab + VaN b + \7 b N a ] 



(6) 



where n is the unit vector field normal to E. One obtains 



S=[dbf d 3 x^ ^(R {3) + K ab K ab - (q ab K ab ) 2 ) + ^L(T - N a d a T) 2 - ^ q ab d a Td b T - NV T (T)+ 



■^{j> - N a d a ct>f - j^da^ ~ NV^) 



(7) 



Knowing that involves only spatial derivatives of q a bi it is easy to compute the conjugate momenta: 



M = ^_ {qabq cd_ q ao qM)Kcd 
oq ab 2k 



5S 
SN a 



= 



SN 



(8) 



PT 



ss_ 

8f 



N 



(T - N a d a T) 



Using these, we can rewrite (0 in the canonical form: 

K ab q ab + PtT + -NC- N a C a 



S= I dt I d 3 x 



(9) 



where 



C 



2k 



Kab^ 



ab 



ab\2 



2k 



2^q' 1 2 
[C a = -2q ac V b ir bc + PT d a T + 7r a a 



+ ^P 2 t + ^q ab d a Td b T + y/qVr(T) + ^ir 2 + ^-q ah d a , ><),,, > - x 7]\'A"'> 



(10) 



Looking at the action (j9|) , one can see that the phase space T and coordinates thereon in which Poisson brackets have 
the canonical form are manifest: we may separate them as 



r = r' r x r T x r 



(11) 



where coordinates (N, N a , q ab , tt, n a , n ab ) parametrize the geometric sector r' G , coordinates (T,pr) parametrize the 
time sector Ft, and (4>,TVd,) parametrize the matter sector Tm- However, as N and N a are non-dynamical, we 
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have four primary constraints (per point x £ E): it — and 7r a = 0. These constraints can be directly solved, 
and the transformations they generate can be gauge-fixed by choosing arbitrary N — N(q a t,,ir ab ,T,pt, 4>,P<j>) and 
N a = N a (q a i,, TT ab , T,pt, <j>,p (/>)'■ the other variables do not depend on them, so they are all gauge-invariant under this 
choice. So the phase space of the system reduces to 

r = r G x r T x r M (12) 

where Tq is parametrized by (q a b, Tt ab ) only. The Poisson structure takes the canonical form in those coordinates: in 
other words, the only nonvanishing brackets are 

{q ab (x) 7 7T cd (y)}=Si c S d b ) S^(x > y), {T{x),p T {y)} = S^(x, y), {#*), n^y)} = S^(x, y) (13) 
Conservation of the primary constraints under the evolution generated by the Hamiltonian 

H = J d 3 x [NC + N a C a ] (14) 

implies four secondary constraints (per point x € E): 

C = 0, C a = (15) 

It can be shown that these constraints are conserved with respect to H, so (C, C Q ) constitutes the whole set of 
constraints. Moreover, their constraint algebra closes, so they form a set of first-class constraints (using Dirac's 
terminology) . 

Observables of the theory are those phase space functions F - called the Dirac observables - that Poisson-commute 
with all the constraints: 

{F, C(x)} = {F, C a {x)} = 0, for all x £ S (16) 
This in particular means that any observable F does not evolve, as it commutes with the Hamiltonian (|14[) : 

j t F={F,H} = (17) 

This so-called problem of time is solved by introducing a suitable automorphism on the Poisson algebra of Dirac 
observables. 

It is possible to show that the constraints (C,C a ) encode an important geometrical feature of the theory: 
diffeomorphism-invariance. Indeed, C a and C generate the action on the phase space of difieomorfisms of E and 
diffeomorphisms off S (in the normal direction) , respectively. For this reason they are often called respectively vector 
constraint and scalar constraint. 



III. BACKGROUND SECTOR OF T 



We want to consider generic linear perturbations on a fixed background. The ideal background (which is also the 
physically meaningful one) presents homogeneous and isotropic space E. We are assuming throughout this paper, 
that E is a 3-torus parametrized by fixed coordinates a; 1 , cc 2 , x 3 € [0, 1) (where the interval [0, 1) is endowed with the 
topology of a circle). Those coordinates will be used to define the "background" sector, a subspace 

r<°> = rg>xrWxrJ°> c r G xr T xr M (is) 

is the homogeneous isotropic part of the geometric sector Tq of the phase space: consists of points (q^ b , nfy) 
such that 

q^(x)=e 2 «6 ab , ^(xH^V (19) 

where a and ir a are constant. Therefore, is freely parametrized by (a,ir a ). 

In r^, the matter field <j> is assumed to be absent, in the sense, that 7r^ = and <fi — (f>o, the minimum of the 
potential V^. For simplicity let us assume that 



0o = Wo) =0 



(20) 
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although in future it may be also be interesting to study the consequences of the spontanues symmetry breaking in 
this context. Therefore, T$ consists of a one point (0, 0). 

Contrary to <j), a nontrivial time field T is necessary. To be consistent with the homogeneity of the space, we choose 
the background T homogeneous as well: this means that the infinitely many degrees of freedom sitting in T(x) and 

Pt{x) are reduced to a unique one. Hence, (T^\p>p) freely parametrize where 

T (0) 



T(x) = TW PT (x) 
The subspace can be intersected with the constraint surface 

r c cr 



(0) 

Pt 



(21) 



(22) 



consisting of solutions to the constraints (fTBj) . Since everything is spatially homogeneous, the vector constraint 
C a (x) = is automatically satisfied for every point of T^°\ In fact, only the homogeneous part of the scalar constraint 
survives (below, we denote the restriction of the constraint C to r' ) by C^): 



C (0) (JV) = / d 3 xN(x)C {0 \x) = e 



3o 



6q Vt(T (0) ) 



d 3 xN(x) 



j ^ '\2 

Therefore, at the intersection with the constraint surface Tc, the points of r' ) additionally satisfy a constraint 
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These points correspond to a FRW spacetime 

g$dx»dx» 



-dt 2 + e 2aW 6 ab dx a dx b 



(23) 



(24) 



(25) 



satisfying Einstein equation with energy-momentum tensor given by the clock field T = T^(t) (which satifies Klein- 
Gordon equation in a spacetime of the form ([25])). The dependence on the variable t of a, n a ,T^ and p^ (upon 
assumption N(t) = 1, and recalling that in the chosen coordinates J s d 3 x = 1) is given by Hamilton equations: 



K 

"6 e " 



3a„ 



*(Pt) 2 



3e ia V T 



(26) 



ji(Q) 



.(0) 

L Pt 



e-^pP 



— 3ct 



dV T 



A solution of this system of equations which also satisfies the constraint (|2"4")l yields the (dynamical) background FRW 
metric gffi on which usual perturbation theory is developed. 

As we are going to see, in our approach to the full theory we do not intersect the homogeneous isotropic sector 
r(°> with the constraint surface Tc, nor we fix the background dynamics to be of the form (|26|) . Instead, all the 
homogeneous isotropic sector r*- ^ will be coupled to the perturbations and together they will obey the dynamics of 
the full theory. Not surprisingly, at linear order in the perturbations, the background dynamics does reduce to (1261) . 

so indeed our homogeneous part of the metric gffj is of the FRW type in that approximation. Nevertheless, it is 
important to notice that: 

(i) In our framework - designed for canonical quantization - the degrees of freedom a,TT a ,T^ and p^? will be 
treated on the same footing as the remaining degrees of freedom; 

(ii) (1241) and (1261) are true only up to the linear order: if one considers higher orders, the back reaction is present, 
and our homogeneous isotropic part will not anymore be a solution to 



IV. COORDINATES ON THE FULL PHASE SPACE 



Now, we go back to the full phase space, and define on it some clever coordinate system, adapted to the background- 
perturbation split (which will be performed in the next section). All formulae appearing in the present section are 
exact and valid for every point 7 € T of the full phase space. 
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A. Extension of the background coordinates to the full phase space F 

The functions parametrizing the homogeneous isotropic subspace T^ ' can be extended to the full phase space. We 
do it in the following way. Given (q ab ,Tr ab ,T,p T7 (j) } p t p) e T, define 



a 






= 2e 2a S ab J s d 3 xir ab 


T (0) 


= k d * xT 


(0) 





(27) 



In this way, a,TT a ,T^ and pip become functions defined on the phase space V of the full theory. Notice that each 
of them is defined globally on T. 3 We call them background coordinates (or background variables) on I 1 . 4 Restricted 
to the subspace T^ ', they coincide with the coordinates introduced in the previous section and denoted in the same 
way. 



B. Perturbation coordinates on F 



Next, for every x e E, on the full phase space T we define functions Sq a b(x), 5n ab (x), 5T(x), Spt(x), 5(j)(x), Sn < p(x) 
which, together with the background coordinates, form a coordinate system on T: 



' 5q ah {x) 


= q a b(x) 


- e 2a 


5ir ab (x) 


= TT ab ( X ) 


ir a 
6 


ST(x) 


= T{x) - 


. T (o) 


6p T (x) 


= Pt(x) 


-pP 


6<p(x) 


= m 




K 8ir ( j > (x) 


= 7T^(x) 





(28) 



Let us call them perturbation coordinates. They satisfy the following identities: 

J d 3 x5 ab Sq ab (x) = J d 3 xS ab 8TT ab (x) =0, J d 3 x8T(x) = J d 3 xS PT {x) = 



(29) 



Such relations constrain the perturbation coordinates. 

Using the background coordinates and the perturbation coordinates, we are simply parametrizing points in T in 
terms of their suitably defined components in 1^°) C T and the rest. It is good to keep the phase space picture in 
mind, since later it will be essential that we do not forget about the background sector (as usually done in standard 
cosmological perturbation theory, by demoting its degrees of freedom to fixed parameters instead of dynamical vari- 
ables). Indeed, while for the classical theory there is no difference, to properly prepare the set up for quantization (of 
the full system: perturbations and background) this is the correct way to go. 



3 Indeed, 8 ab q a i,(x) > ensures that even a is defined globally. 

4 This terminology comes simply from the fact that, in the next section, we will be Taylor-expanding certain phase space functions 
(namely, the constraints) around a ncighbourghood of points (a, 7T a , T(°), 0, 0) £ T. As already explained, no dynamical propery is 
taken into account here. 
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C. Fourier mode-decomposition 



The coordinates (x a ) = (x , x 2 , x 3 ) fixed on S can be also used to introduce a mode-decomposition (or real Fourier 
transform) of fields. First, let us define the usual Fourier transform (with respect to coordinates (x a ) for S) of a field 
f(x) by 



/(fc) = / d*xe-* k x f(x) 



(30) 



We think of A; as a spatial vector (i.e., tangent to £), which labels the mode of /. It takes values in the lattice 
C = (2ttZ) 3 . 

For a real- valued / as our fields, it is f(—k) — /(fc)- Thus, to isolate the truly independent modes we work with 
the real Fourier transform. Split the lattice C into "positive", "negative" and "zero" vectors: 



C+ = {k e C : (fci > 0) V (fci = A k 2 > 0) V (h = k 2 = A fc 3 > 0)} 
£- = {k e C : (fci < 0) V (fci = A k 2 < 0) V (fci = k 2 = A fc 3 < 0)} 
Co ={0} 



(31) 



Clearly, we have C = £+ U £_ U Cq. Then, we define the real Fourier transform of f(x) as 

/(*0 : 



(f(k) + f(-k))/V2 iffce£+ 
(f(k) - f(-k))/iV2 if fc g £_ 
/(0) if fc = 



(32) 



More explicitely, we have 



V2 



J d 3 x (e ik - x + e~ lk - x ) f(x) like 



/(*) = 



/ d 3 x (e lk - x - e- lk ' x ) f(x) 



V2 

{ fd 3 xf(x) 



if fc e C- 
if k = 



(33) 



Knowing /(fc), we can easily reconstruct the field: from Fourier antitransform, we write 

/(x) = /(0)+ e ik ' X f(k)+ e* x /W = /(0)+ £ [e* a /(*) + e~ ik - x f(-k) 



(34) 



fee £4 



fce£_ 



where we used the fact that fc € £- is quivalent to — fc G £ + , and then replaced the dummy index — fc with fc. Inverting 
the definition (f3"2"j) of /(fc), one has 



V2/(fc) = /(fc) + /(-fc), iV2f(-k) = /(fc) - /(-fc), for fc e £+ 
thus, by adding and subtracting these two, one finds respectively 

f(k) = -^(f(k)+if(-k)), f(-k) = ^=(f(k)-if(-k)), forfce£ + 



V2 



V2 



Replacing these in (|34l) . we finally obtain the inverse transform: 



m - /(o) + 



fee £4 



V2 



*/(-*:) 



V2 



(35) 



(36) 



(37) 
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Specifically, for our fundamental fields the mode-decomposition is the following: 5 



q ab (x) =e 2a S ab + 5q ab (0) + Y: 



5q a b(k) 



V2 



+ i5q ab (-k) 



V2 



ir ab (x) 



(i 



-s ab + Sir ab (o) + J2k 



keCj 



SlT ab (k) 



V2 



+ i57T ab (-k) 



V2 



T(x) =T(0)+j: keC 

p T {x) =Pt(0) +Efc S £ + 

<t>(x) = 54>(o) + £ fce£+ 



ST(k) 



V2 



+ iST(~k) 



V2 



6p T (k) 
6<j>{k) 



V2 - + V2 



V2 



+ i<y^(-fc) 



72 



(fc) 



V2 



+ iSn^—k) 



V2 



The fc = mode corresponds to the homogeneous part. So we have 

T(0) = T(°), M0)=*4 0) 
For the metric perturbations, the /c = case is nonzero only as a traceless matrix: 

S ab Sq ab (0) = 0, <5 a6 <57f Qb (0) = 



(38) 



(39) 



(40) 



As a confirmation of this, notice that the k = mode is by definition /(0) = / d 3 xf(x). But since (|29|) holds, we 
have directly the constraints just mentioned. 



D. The scalar, vector and tensor modes of the metric 

The treatment of the metric and its momentum requires some more work. Indeed, since 8q ab {k) defines a symmetric 
3x3 matrix per each k, once we fix the mode k we can expand Sq ab (k) on a basis for the 6-dimensional vector space of 
simmctric 3x3 matrices. A good basis in the space of symmetric bi-covariant tensors is {^4^,} (m — 1, 6) defined 
by 



Kb = Sab 



A 2 

A ab 



k a k b 1 
~fc^ ~ 3 ab 



A lb = ~^ kaVb + kbVa ^ Kb 



Kb = -^{v a w b + v b w a ), A 6 ab 



1 

V2 



(k a Wb + k b W a ) 
(v a v b - W a W b ) 



(41) 



where v and w are spatial vectors forming with k an orthogonal basis of the momentum space R 3 (with respect to 
the fiducial metric 5 ab , which is also used to raise and lower indices for v, w and k). The normalization of v and k is 
chosen to be v 2 = w 2 = 1/k 2 . 

The subspaces spanned respectively by (A 1 , A 2 ), (A 3 , A 4 ) and (A 5 , A 6 ) are said to comprise the scalar modes, the 
vector modes and the tensor modes. They have the following properties: 



5 Since we are using the Fourier transform with respect to the fiducial metric, <5 a (,, the wavevector k a is not the physical momentum. To 
get the physical momentum, one needs to mutiply by the inverse of the scale factor. 
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• tensor matrices satisfy k a A" l b (k) = 

• vector matrices satisfy fc a fc b A^(fc) = 

Also, notice that all matrices except from A\ h satisfy 

5 ab A™ b (k) = (42) 

We decompose 8q ab (k) in this basis, 

Sq ab (k) = q m (k)A2(k) (43) 

where q m (k) denotes the mth component, for m = 1, ...,6. 
Similarly, one can expand 5TT ab (k) on the dual basis {i^ 6 }: 

8% ab {k)=p m {k)A ab {k) (44) 



The dual basis is given by 



1 3 f h a h h 1 

Af = -5 a \ Af = -(^ r - -6 ab 

1 3 ' 2 2 V k 2 3 

Af = — (k a v b + k b v a ), Af = -j=(k a w b + k b w a ) (45) 



k 2 k 2 
Af = -=(v a w b + v b w a ), Af = -=(v a v b - w a w b ) 
v 2 v 2 



It is easy to check the duality, i.e. that Tr(A m A n ) = AfA^ a = AfA n ah = 5™ . Moreover, all these A matrices are 
normalized with respect to the scalar product induced by the fiducial meric 5 a b, namely 

(A, A') = 5 ac 5 bd A ab A' cd , (A, A') = S ac 5 bd A ab A' cd (46) 

except for A 1 , A 2 , A±, A2, for which we have 

{A\A X )=Z, (A 2 ,A 2 ) = 2/Z, (A 1> Ai) = l/3, (A 2 ,A 2 )=3/2 (47) 

These scalar matrices are left non- normalized to keep the agreement with formule in [14j . Note however that in our 
case the matrices A do not involve any dynamical variable (in particular, they do not depend on a), so no nontrivial 
Poisson algebra is hidden in expansions (|43|) and (|44j) : everything is contained in our new variables q m (k) and p m (k). 
Finally, note that it is always possible to choose v and w in such a way that A are symmetric under k — > — k: we do 
this, so in the following we will have 

KM = AZ(-k), Af(k) = Af(-k) (48) 



E. Summary of the section: resulting coordinates on F 

Let us summarize. We have introduced on the full phase space V the following system of coordinates: 

• the background coordinates: 4 numbers (a, ir a , T(°\pj! ), 

• the perturbation coordinates (homogeneous part, that is, for k = 0): 12 numbers (Sq a b(Q), <57r ab (0), 8(f)(0), 5^(0)), 

• the perturbation coordinates (inhomogeneous part, that is, for k 6 C — {0}): 16 numbers per each fc, 

(q m (k),p m (k),Sf(k),6pT(k), 6$(k), Sn^k)). 

Again, we remark that this terminology is in function of the next section (and the remainder of the paper), but up to 
now the "perturbation components" are defined for every point in phase space T, and are finite. 

The fundamental Poisson algebra of the original ADM variables straightforwardly induces the following Poisson 
algebra on the new ones: 

{a,7r Q } = l, {T(°\pW} = l, {Sq ab (0),Sn cd (0)} = 895* - \5 cd 5 ab , {^(0), 5^(0)} = 1 

(49) 

{q m (k),p n (k')} = SlS kX , {6f(k),Sp T (k')} = 5 k ,v, W(fc),,57r (fc')} = Sk,k> 
for fc, k' ^ 0, and the remaining Poisson brackets vanish. 
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V. THE CONSTRAINTS UP TO THE FIRST ORDER 



A. The expansion 



We now turn to the constraints of the theory. At this point, it is convenient to expand them for "small" perturbation 
variables (|28|) . That is, we are applying the Taylor expansion formally given by 



+ Sf) = F(/<°>) + F'(fW)6f + l -F"(f^)6f 2 + 0(6 f) 



(50) 



In our case, / stands for the fields, for the background variables, and 6f for the perturbation variables. The 
decomposition / = /(°) + Sf reads 



r q a b{x) 


= e 2a 5 ab + Sq ab (x) 


7T ab (x) 


6 


T(x) 


= T(°) + 8T(x) 


p T (x) 


= pf' + Sp T (x) 


4>{x) 


= 5<f>(x) 


< n<j>(x) 


= <57t,j (x) 



(51) 



In light of the decomposition (1511) . we are to expand accordingly the scalar and vector constraints: 

C(x) = C<°> (x) + C« (x) + (x) + 0(6 f), C a (x) = tf<°> ( x ) + C« (x) + Cf> (x) + 0(6f 3 ) (52) 

where (and C^) collects all the terms which are 0th order in the perturbation variables 8f (i.e., the first term 
in the Taylor expansion), (and C^) collects the 1st order terms (the second term in the Taylor expansion), and 
(and C^) collects the 2nd order terms (third term in the Taylor expansion). The idea is thus that we retain the 
full constraints and deal with their exact solutions, but we write such solutions explicitly only up to the first order in 
the perturbation variables. More precisely, we have that the first order expansion of the solutions satisfies 

C^(x)+C^(x)=0(Sf 2 ), Ci \x) + Ci x \x) = O(Sf) (53) 

Therefore, in this linear approximation we will solve 

C {0) (x) + C {1) (x) = 0, (x) + C { a 1] (x) = (54) 

and ignore C(") and C ( a n) for n = 2,3,.... On the other hand, we will retain to describe the dynamics at first 
order. 

Now, using (|33|) . we write the real Fourier transforms of the constraints: 



C(0) 
C(k) 

C(k) 

C a (0) 

C a (k) 
C a (k) 



fd 3 x + fd 3 xCM(x) + fd 3 xCW(x) + 0(6 f 3 ) 

1 



/ d 3 x (e lk ' x + e~ ik - x ) [C(°) + C« (x) + (x)] 



-j= f d 3 x (e lk - x - e~ ik - x ) [C(°) + C« (x) + (x)] 
Ci 0) fd 3 x + f d 3 xC a [ 1) (x) + f d 3 xd 2) (x) + 0(6 f 3 ] 



0(6f 3 ) 
0(6f 3 ) 



if k e £+ 



if k e C- 



(55) 



V2 



V2 



f d 3 x(e lk - x +e~ %k - x ) C { a ] +C ( a ) (x) + C ( a ) (x) + 0(5f 3 ) iffce£+ 



.(1) 



(2), 



Id 3 



) Ci 0) + Ck 1, (x)+C^(x) 



(2), 



0(6f 3 ) if k e C. 



These contain many terms, but the following remarks will help us simplifying them: 
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• As already pointed out above, vanishes identically. 

• We already said that we disregard the second order terms except for the dynamics. Since as far as H is concerned 
we are free to choose the lapse and the shift, we will select N(x) — 1 and N a (x) = in (fT4|) . thereby obtaining 
simply H = J d 3 xC(x) = (7(0). It follows that we must retain the second order only in (7(0). 

The terms fd 3 xC il) (x) and $&xC£\x) can be seen to vanish because of (|29[) . Also, one sees that 



C^fd 3 x(t 



±e 



— ik-: 



C^Sk,o, so it vanishes for all k ^ 0. 



Applying these remarks, we see that: 

(i) C a (0) — 0(6 f 2 ) is identically satisfied 

(ii) using (|23|) . the constraint (7(0) reads 

(7(0) = 



-3ct 



0(6f 2 ), 



which will later be used to solve for function of the other background variables. 



(iii) we are left with 4 constraints per each k 7^ 0: 



C(k) 
C(k) 
C a (k) 

C a (k) 



1 



-^fd?x(e ik - x + e- ik - x )CW(x) + 0{6f 2 ) 



V2 



■ x )C<- 1 \x)+ 0(Sf 2 ) 



/ d 3 x (e ik ' x + e~ lk - x ) C ( a\x) + 0(5 f 2 ) 

4= / d 3 x (e ik - x - e~ ik - x ) C ( a\x) + 0(5 f 2 ) 
V2 



if k e £+ 
iikeL- 
if k G £+ 
iikeL- 



(56) 



(57) 



B. Explicit form 

In order to find their explicit form, we need to first compute the linearized constaints as function of a space point x, 
C^(x) and Ca(x). Plugging the decompositions ([5T|) into (fTUl) and keeping only the terms linear in the perturbation 
variables Sf, we find 6 



-3a 



P ( tSpt + < —V T (TW)q$ ) 6q ab + e 6a V+(T^)5T 



4 V 18 



. (>0 



(P 



(0)x2 \ „ab£„ t K7F " _(0) rab 



(58) 



{ C ( a L> = ir^OJqtc - 2q™d c 6n bc - 2n b ( ^d c 5q ab + p^dJT 

In obtaining these linearized constraints, we used some nontrivial facts: 

• The determinant q gets contributions only from the diagonal terms (since at 0th order the non-diagonal terms 
are zero), so one has 



<l 



(o) 



d 9ll<722 <?33 + 111 0922933 + 9n 922 "1 



JO) 



(°) _l „(°U°); 



33 



^6a 



+ e 4a <5 ab ^ afc 



(59) 



We dropped all the terms coming from the matter field 0, since they always involve a 0th order factor, i.e. either </>(°) or 7r^ ' , which 
are zero by definition. In this sense, the role of of a "test field" is mathematically justified. 
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For the spatial Ricci scalar we used the relation 

= q °» ( dc T a c b - d b r a c c + r a \r c c d - r a d c r b c d ) 



(60) 



However, the last two terms are second order in Sf because they involve products of two 1st order objects 
(namely, the spatial derivatives of q ab sitting in Christoffel symbols, which are 1st order because the 0th order 
metric is homogeneous). So one is left with 



R {3) = (<?(o C )9(o) - 9(oV(o)) d a d b Sq cd + 0(5 f) 



(61) 



• Since n is a tensor density of weight 1, its covariant derivative is 



V^bc o be i -p b dc i -p c bd p rt ^ 
Q 7T = (J a TT +l ad 7T +l ad 7T ~ 1 d a 7T 



d „bc 



(62) 



Notice that in Vb7r which appears in the vector constraint, the last term cancels with the second one. 

At this point, we can compute the real Fourier transforms (1571) . To sec how it works, consider the first one, C(k) 
for k € £+. It is: 

E(k) := C {1) (k) = l=Jd 3 x (e lk - x + e~ lk - x ) C {1 \x) = 



-3a 



V2 



-3a 



V2 



^ - d a d b Sq cd -\(^ + ipP?) <$)S gab 



-q a >5n ab +p^5 PT + C —V T (T^)q^Sq ab + e e «V'(T^)6T 



ik-x | ik-x\ 



1 / KTTi 



4 V 18 



-^i°> ab +P [ t ] Spt + ^V T (T^)q^6q ab + e e ^(T^)ST 



-3a 



h' 7V, 



^ («(oV(0) - 9(0^(0)) kak b Sq cd + (P^) 2 ) 



. (in 



-f^S^+p^Spr + C —V T (T^)q^6q ab + e to t£(T<°>)*f 
where in the last step we used (|33[) with respect to the perturbation variables. We can rewrite this as 



(63) 



E(k) 



2 k 



(k 2 6 ab -k a k b ) AZ(k)q m (k) 



- 5 a. 
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f + (p^) 2 - 2e 6a V T (T^) ) S ab AZ(k)q m (k) 



KTT a e 



-S ab A^(k)p m (k) + e' 3a p { ° ] 5pT(k) + e 3a V^(T^)Sf(k) 



(64) 



having replaced the explicit expression for the background variables, and having expanded the perturbation variables 
of the metric (and conjugate momentum) in the {.A^} basis. Repeating the computation for C(k) in the k £ £_ 
case, we obtain the same object, which is then regarded as the scalar constraint satisfied by each mode k E £ — {0}: 



E{k) = - ^— ( ^ + [pPf - 2e^V T {T^)\ qi (k) - —k 2 qi {k) + e —k*q 2 {k) 



18 



3k 



Kir a e 



l -p l {k) + e- ia p { p8p T {k) + e 3a V+{T^)5f(k) 



(65) 



where we used (l4"TT) and P5|) to write E(k) as a function of the dynamical variables only. 

One obtains C^\k) in a similar way. Moreover, since C a (k) really encodes 3 constraints, we need to separate them. 
To do this, we project C a (k) along the 3 orthogonal vectors k, v and w at our disposal: 



M(k) = k a C a (k), V(k) = v a C a (k), W(k) = w a C a (k) 



(66) 
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Explicitely, these 3 constraints satisfied by each mode k € C — {0} are given by 

' M(k) = ^f^gi(fc) - ^^q 2 (k) - ^V(A) - 2e 2 V(fc) +p(?^f (fc) 
6 9 3 

< V (k) = ^£^ 93 (fc)+2e 2 V(fc) (67) 

W(jfc) = H°l 2a g 4 (fc)+2e 2 > 4 (fc) 

Notice that, apart from the background variables, the linearized constraints E'(fc) and M(fc) only involve the scalar 
modes (namely qi(k), q^ik), p 1 {k) 1 p 2 (k), ST(k) and 6pr(k)), whereas V(k) and W(k) only involve vector modes 
(namely qa(k), qi(k), p 3 (k) and p A (k)). It follows that the linearized constraints E and M only constrain the scalar 
sector, while V and W constrain the vector sector. Interestingly, the tensor sector is left completely unconstrained. 



C. Preliminary analysis of independent degrees of freedom 

Recalling that each constraint reduces the number of degrees of freedom by 2 (1 for the reduction onto the constraint 
surface, and 1 for fixing a gauge - or equivalently for identifying each 1-dimensional orbit with a single point), we can 
then proceed with the counting of degrees of freedom: 

• For k = 0, up to the first order the constraints (7(0) and C a (0) constrain only the background coordinates 
a,TT a ,T^ and p§ by ((56)) . That constraint can be solved for p^' and used to gauge fix T^°K On the other 
hand the traceless variables 5q a b(0), i57f a6 (0), and the variables </>(0), 7^(0) are unconstrained. 

• For every k ^ and given background coordinates a, 7r Q , and pf \ the scalar sector of the phase space is co- 
ordinatized by (qi(k),p 1 (k),q2(k),p 2 (k),ST(k),SpT(k), S^>(k),6Tt^(k)), so it has dimension 8. On it there are the 
2 constraints E(k) and M(k), so the dimension is reduced by 2-2 = 4. We conclude that the corresponding sector 
of the reduced phase space has dimension 8 — 4 = 4: in other words, there are 4 gauge-invariant (i.e., physical) 
scalar degrees of freedom. 2 of them can be chosen to be (S^(k),dn^,(k)). The other 2 independent degrees of free- 
dom can be chosen to be particular functions of the remaining variables (qi(k) , p 1 (k) , q2(k) , p 2 (k) , ST '(k) , SpT(k)) 
(see later). 

• For every k ^ and given a, 7r Q , the vector sector of the phase space is coordinatized by 
(q3(k),p 3 (k),q4(k),p 4 (k)), so it has dimension 4. Imposed on those variables there are 2 constraints, V(k) 
and W(k), so the dimension is reduced by 2 • 2 = 4. We conclude that the reduced phase space has dimension 
4 — 4 = 0: in other words, vector modes are completely non-physical, and can be gauged away. 

• For every k ^ and given background coordinates a, 7r Q ,T(°) and p?\ the tensor sector of phase space is 
coordinatized by (q§(k),p 5 (k), gg(fc),p 6 (fc)), so it has dimension 4. Up to the first order, there are no constraints 
imposed on those variables, so there is no reduction in dimension: we conclude that the reduced phase space 
has dimension 4, i.e. there are 4 degrees of freedom. Obviously, they are (q^(k),p 5 (k),qQ(k),p 6 (k)) themselves. 
These (or rather, the 2 configuration variables q$ and q 6 ) are the 2 polarizations of the graviton. 

For the reader familar with the perturbative approach to the canonical gravity often used in cosmology, it is important 
to notice the difference of our approach. In the standard approach, one often introduces functions of the perturbation 
coordinates which Poisson-commute with the linearised constraints E, M, V, and W. In that approach they are 
Dirac observables and play a fundamental role. In our approach, the transformations generated by the 2nd order 
terms of the Taylor expansion of the full constraints also contribute to the gauge transformations. Moreover, for 
us also C' 2 )(0) is a generator of gauge transformations. Therefore, Poisson commuting with E, M, V, and W is 
not sufficient to be gauge invariant in the sense of the current paper. The consequence is that, within the approach 
presented in this work, gauge-invariant observables (the Dirac observables) are not just simple combinations of the 
perturbation coordinates on the phase space even up to the first order. Having said this, in the following we may 
use the terminology "gauge- invariant variables" for those coordinates that commute with E, M, V, and W. One 
should however bear in mind that they are not gauge-invariant in the full theory, and in fact we will be carrying out 
gauge-fixing in order for them to represent observable quantities. 
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VI. THE CONSTRAINTS UP TO THE FIRST ORDER: SOLUTIONS AND GAUGE-FIXING 

A. Solution to the Constraints 



In this subsection, we solve the constraints up to the first order. Every solution to the full constraint coincides 
with one of our solutions up to the second (or higher) order in the perturbation variables. Specifically, we will show 
the general solution to all the constraints, thus reducing to the constraint surface Tc C T. Next, we will choose a 
family of slices of Tc transversal to the orbits of the gauge group generated by the constraints. We will finally find 
field variables freely parametrizing the slices and use them to fix the gauge, thus reducing to the physical phase space 

Tphys • 

To strat with, the constraint (|56p up to the first order reads 



Wt) 2 + e 6a V T (T^) 
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We solve it with respect to Pj! > finding 



(o) 
P T 



± 



^wl-2e**V T (T«») 



The gauge transformations generated by C (0) can be used to fix arbitrarily the value of T^°> : 

T (o) _ T = o 



(68) 



(69) 



(70) 



The vector constraint equations 

' V(k) = 
W{k) = 

can be immediately solved for the momenta p 3 and p 4 : 



for aU k € C - {0} 



p 3 (k) 



-4a 



6 



-q 3 (k), p\k) 



ir a e 



-4a 



6 



-q 4 (k) 



(71) 



(72) 



The associated configuration variables, ((73 (fc), qi(k)), are free. They can be used to parametrize the gauge-orbits 
of the constraints v a C a (k) and w a C a (k) (which, at linear order, are nothing but V(k) and W(k)) in the constraint 
surface Tc C T. 

Let us now consider the scalar constraint equations: 



E(k) = 



for alike C- {0} 



(73) 



M(k) =0 

We proceed as above: we simply solve them for the momenta p 1 and p 2 . The result is the following: 

3e 4a 



p\k) 



-V^ST(k) + Sp T (k)- 



Tr a e~ 4a 3fc 2 9(p { ° y ) 2 e- 4a 9e 



2d 



2kit, 



-V t 



qi(k) 



V 2 {k) 



-14 



ST(k) - ?Z Sp T (k)- 



(74) 



7r„e 



-4<:> 



3(^ 0) ) 2 e- 4 « 

4k7T q 



2d 



3e 

2 KIT, 



-V T 



qi(k) 



ir a e 



-4a 



3k 2 7T q 



92 (A) 



1G 



In this way, (qi(k),q2(k),5T(k),6pT(k)) can be used as free coordinates on the constraint surface Tc- Two of these 
free variables should be chosen to parametrize the gauge-orbits of C(k) and k a C a (k). The remaining two functions 
will represent the "physical" degrees of freedom. 

Finally, as for the tensor sector, we notice that it is completely unconstrained. Therefore, all variables 
{Q5{k),P 5 {k),qe(k),p 6 (k)) are free: they are the "physical" degres of freedom associated with the graviton. 
The remaining variables (namely, the traceless matrices Sq a b(0) and c)7r afc (0), as well as all the modes 8(f>{k) and 
8ir^(k)) are also free. 

Of course all the "free" functions listed above are still subject to the gauge transformations, hence they represent 
the physical degrees of freedom in a gauge-dependent way. 

B. Gauge-Fixing 

At this point, we are following the so-called Reduced Phase Space formalism, in which one first solves the constraints, 
and then identifies each gauge-orbit with a point in the physical phase space. This procedure is generally regarded 
as ideal, but can be implemented explicitly only in few cases. Indeed, it is usually a hard task to identify a set of 
gauge-invariant functions to coordinatize the physical phase space, and moreover they have usually a very complicated 
Poisson algebra, which makes canonical quantization practically impossible. 

A way out of this problem is the so-called gauge-fixing procedure. One chooses a set of functions to play the role of 
gauge-parameters, i.e. parametrizing the gauge-orbits of the constraints. Points along the same orbit are physically 
equivalent, so we can simply choose one to represent that specific orbit (physical state of the system). This amounts 
to fix the value of the gauge-parameters, choosing a "slice" in the constraint surface which mirrors the physical phase 
space. This slice is endowed with a symplectic form by simply pulling back the kinematical symplectic form along the 
embedding, and hence it is a good phase space. It is coordinatized by the remaining variables, which thus represent 
the physical degrees of freedom of the system. More precisely, the interpretation of these surviving variables is the 
following. We imagine to be dealing with Dirac-observables (i.e., functions that commute with all the constraints), 
denoted 7 ,g- Now, j 1 is the value of C 7j g when restricting to the chosen slice (here collectively denoted by the 
gauge-fixing conditions G a — 0). There is a well developed formalism (2l| - [27| to treat these objects, and we will be 
using part of it when studying the dynamics. For now, it is enough to know that gauge-fixing is a "legal" procedure 
to reduce to the physical phase space. 

In our case, we have to deal with 4 constraints C(k), C a {k) for each k ^ and with C*(0). Good gauge-parameters 
for such constraints can be chosen to be (qi(k), q2(k), qs(k), q\ (k), T^). The gauge-fixing conditions that we choose 
are the following: 

qi (k) = q 2 (k) = q 3 (k) = qi(k) = 0, for all k £ C - {0} (75) 

T(°) - r = 0, with tg! (76) 

Note that the value r of gauge-parameter is left free: it will be used to label the Dirac observables and its changes 
will be used to describe their dynamics. 

VII. THE PHYSICAL PHASE SPACE, OBSERVABLES AND THEIR DYNAMICS 

A. Physical phase space and observables 

Mathematically, the physical phase space r p hy S for the theory we are considering in this paper is the space of the 
orbits in the constraint surface Tc of the gauge transformation group generated by the scalar constraints C (x) and 
the vector constraints C a {x). The space r p h ys can be embedded in Tq as a slice which intersects each orbit exactly 
once. This is exactly what we did up to the first order in the previous section by solving the constraints and fixing a 
gauge. Denote the image of the embedding by Fp hys , that is 

r phys -> r^ hys c r c c r (77) 

where r is the parameter used in the gauge conditions. The surface Fp hys can be parametrised by the following 
variables originally defined in all the kinematical phase space T, 

( 7 j) = (a,n a ,5q ab (0),5t ab (0), 50(0), 5^(0), 6f(k),8p T ik),q 5 {k),p 5 (k),q 6 {k),p e (k),S4>(k),5^(k)) , (78) 

for all k £ C - {0} 



17 



The embedding (|77|) determines the surface r p h ys in T up to the 1st order as follows: 

Ql (k) = q 2 {k) = q 3 (k) = q 4 (k)=p 3 (k)=p i (k) = 



pi(fc) = (p^ 0) <5pT(fc) + e ea V^6T(k) ) , P 2 (k) = -Et^ 



(o) 

"".J--." "■.i-'A-TY/.i'i - - -^J ,5p T (fc) + e- 2Q \ Vj y-—Vt ) <5T(fc) 



6 a 



e 

K7T, 



(79) 



T(°) = r, ?4 0) = i^TT^ - 2e«°y T (T(0)) 

The pullback of the coordinates (|78p to r p h ys defines coordinates 

( 7 J) = (a T , <, 5q ab (0y, 6n ab (0Y,Sm T , «r*(0) r , 5T{k)\5p T (k)\ q 5 (ky,p 5 (ky, q 6 (k) T , p 6 (k) T , 84>(k)\ Sn^ky) , 

(80) 

for all fc € £ - {0} 

on r p hys- Note that, as emphasised by our notation, each of the functions 7J depends on the fixed value of r. We 
will come back to that dependence below, while defining the dynamics. 

What is independent of the embedding ((77)) is the symplectic form fi p hy S , i.e., the pullback to r p h ys of the symplectic 
form in T. This is the physical symplectic form. To explicitely find it, we simply pull back the kinematical symplectic 
form 

n = da A dn a + eff (0) A dp { ° } + ^dSq ab (0) A d<57r afc (0) + dS^>(0) A d<57r (O)+ 



E 

fce£-{o} 



d5f(k) A dSp T (k) + d 1m(k) A dp m (k) + d54>(k) A dSn^k) 



(81) 



Owing to the convinient choice of gauge-fixing conditions, the pullbacks of dT^ , dq\ , dq 4 by (|77[) vanish identically. 
Therefore the pullback of the symplectic form reads simply 



Q phys = da T A d-K T a + -dSq ab (0) T A dSn ab (0) T + d6^(0) T A d^(0) T - 



E 

kec-{o} 



dST{k) T A d5p T {k) T + d< lra{ky A dp m (k) T + dS4>(k) T A dSn^ky 



(82) 



The Poisson algebra that r2 p h ys defines on r p h y s is easily found: 

{« T ,<W = 1. {Sqab(0) T ,Sn cd {0y} phys = Sffi - \5 cd 5 ab , {^(0) T , <57f (O) T } phys 

{q 5 (ky, P 5 {k'y} phys = s kt v, {q 6 (ky, P 6 {k'y} P hys = s k , k , 
{sf(ky, sp T (k'y} phys = 5 k<k < , {S4>(ky, 6n^k'y} phys = s ktk , 



(83) 



This is the canonical Poisson algebra, the simplest we could hope to obtain: encouraging fact, in light of the future 
canonical quantization. 

Thus, we conclude that the reduced (physical) phase space r p h ys is coordinatized by the functions (l80l) . Each of 
the variables (|8T)]) defines a Dirac observable. Conversely, every Dirac observable can be represented by a function 
/(7J). The physical Poisson bracket between two such observales can be calculated from (I55|) . This concludes the 
characterization of the kinematical structure of the physical degrees of freedom of the theory. 



B. Dynamics in r p h ys 

The dynamics of the theory is encoded in the dependence of the variables (|80[) parametrising r p h y s on the gauge 
parameter r. Since the Poisson algebra (1831) is canonical for every r, the dependence of the variables on r is a flow of 
canonical transformations generated by some T-dependent function ft phys defined on T p h y s, and such that 



d 

7k 



ii = {7]>ph yB W (84) 
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We call this function a physical Hamiltonian [28|, [29| . Obviously, it does not have the form of the canonical Hamilto- 
nian J d 3 xN{x)C{x) + N a (x)C a (x), because the canonical Hamiltonian vanishes identically on Tq in which r p h ys is 
embedded. On the other hand, h p h ys must follow somehow from the canonical dynamics J d 3 xN(x)C(x)+N a (x)C a {x). 
Therefore, to derive the physical Hamiltonian we go back to the constraint surface Tc in the kinematical phase space 
r, and even to T itself (because Tc is not equipped with a symplectic form or with the Poisson bracket). 
Let us use the projection 

n : r c -> r phys (85) 

to pullback every function / defined on r p h ys , to a function on Tc, 

O f = n*/ (86) 

called Dirac observable. It is constant on each orbit of the gauge transformations, and extended arbitrarily to the 
kinematical phase space T. This Dirac observable Of weakly Poisson-commutes with the constraints, so in particular 

Of, J d 3 xN(x)C(x) +N a {x)C a {x)^ =0 (87) 

Conversely, every Dirac observable defines a function / on P p hys- The kinematical Poisson algebra of the Dirac 
observables corresponding to the variables parametrising the physical phase space Tphys is consistent with the physical 
Poisson algebra in r p hy S) that is 



{O f ,Of,}\r c = { /,/, W |r o 

At this point, let us choose a trivial potential Vt for the theory (|3]): 

V T = (89) 

It follows that the theory is invariant with respect to the translation 

T^T + t (90) 

Therefore, if O is a Dirac observable, so is the following function O t : 

O r {T,...) := 0{T-t,...) (91) 

where stands for the remaining variables which are the same on the both sides of the equality. In particular, 

given 7J, and 7j +Ar , the corresponding Dirac observables C 7 j and Q^t+&t are related as follows, 

7;+ Ar(T,...) = 7j (T-Ar,...) (92) 

It follows that 

d d 

-r<V = — =-7sr<V (93) 

But the right hand side can be calculated from equation (|87|) . where we have still freedom in choosing the lapse 
function and the shift vector. To find the most convenient choice, let us write 

C(x) = -±={ PT {xf - h{xf) = PT [ x) l^ X \ PT (x) - h{x)) (94) 

(we can always do so, since C(x) — and p^/2y/q > imply that C(x) — p^/2^/q < 0, and can then be written as 
the negative of a square). The function h(x) 2 can be explicitely calculated from the scalar constraint, and it has the 
form 

h 2 = -4k (n ab ir ab - ±(q ab n ab ) 2 ^j + -qR^ _ ^ _ qq ab d^d^ - 2qV <p (cf>) - qq ab d a Td b T (95) 
It is immediate to see that h 2 (and hence its square root h) does not involve nor p^: 



Now, a good choice for lapse and shift is 7 

N = 2^/( PT + h), N a = 
In this way, the canonical Hamiltonian H present in equation (|87p. reduces to 

H = J d 3 x(p T (x) - h{x)) = p { ° } - h 
having used the fact that J d 3 xpr{x) = pj! and having defined 



h := I d xh(x) = n*ft, p h ys 



From all this, it follows that 



Note that pi? is a Dirac observable. But (|98|) implies that so is h, and hence 

h = O hphys 

where /i p hy S is a function defined on T^w by h. Finally, we have 



= {O^JJp-h} = 



{0- (J ,h} 



dr 11 



i^phys}phys 
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(97) 



(98) 



(99) 



(100) 



(101) 



(102) 



Purpose of the next subsection is to explicitely express the physical Hamiltonian as a function of the free coordinates 
(7J) on the physical phase space: 



lphys(7/) = h (jJ,Qn = 0,T (0) = T, P n = p n (jj),p? = P^dl)) 



(103) 



where p n (7j), ... (with n = 1,2,3,4) is given by dropping the subscripts r, using (1791 and restoring the subscripts r 
again. 



C. Explicit form of the physical Hamiltonian 

The derivation in the previous section - a self-contained construction which is a special case of the powerfull theory 
of relational observables piTt^ - is exact. However, to explicitely express (|103j) . we need to go back to the expansion 
in the perturbation variables. 

Let us consider the argument (|95j) of the square root. Plugging the expansions (|5Tj) in it, we find 

h 2 = -in [(e 2a S ac + 5q ac )(e 2a S bd + 5q bd ) {^-e~ 2a S ab + Sn ab ) (^- e - 2a 5 cd + Sn cd ) ~ 

_ I ( {e ^ 5ab + Sqab ) ( ^ e - 2a 5 ab + 6n ab )) 2 } + -{S^q + S^q + 5^q)(S^R^ + S^R^ + jWflP))- 
2 V V 6 / / k 



SiA - (S^q + 6<Uq + S^q) \(5^q ah + S^q ab + 5^q ab ) (d a 4>d b <l> + d a T8 b T) + 2V^) 



(104) 



Here, we denoted by S^q, and S^q ab the ith order of q, R^ and q ab respectively. The expansions of such 

phase space functions can be found in the appendix. 



7 Another choice for the lapse function would be N = 2^/q/(px — h). In this case, H = p T + h. But since h > 0, from H = it would 
follow that ' < 0. This is mathematically acceptable (it corresponds to a contracting FRW universe, but is not the physical universe 
in which we live (which is expanding, hence > 0). 
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Expanding these products and keeping up to 2nd order, one finally groups the various terms according to their 
order. Formally, h 2 is of the form 



h 2 = A + Be + Ce' 

where e and e' are respectively linear and quadratic in the perturbation variables. Explicitely, 



A = 



6 



(105) 



(106) 



Be 



S ab Sq ab + 



2K?r„e 2a 



, 5 ab 5n ab + d a d"Sq ab S ab d e d e 5q ab 

6 K K 



(107) 



Ce' = -- 



1 



1 



S ab S cd _ _ S ac s Mj _ ^4* _ j Sl^S***- 

2k7t Iktt 1 1 

— -5q ab 8v ab + —^S ab 8 cd 5q ab 5n cd + -5 ab Sq ab d c d d 6q cd - -5 ab 6q ab 6 cd d e d e 5q cd - 

6 6 K K 

-S ab 5q bc d c d d 5q da + -5 ac S bd Sq ab d e d e Sq cd + -Sq ab 5 cd d a d b 5q cd + -6 cd d a Sq ab d b 5q cd - 

Kj Ki hi k 

-Ls ab 8 cd d e 5q ab d e Sq cd - ±-S ab d d 5q ac d c Sq bd - U cd d a 8q ac d b 5q bd + l-5 ac 6 bd d e Sq ab d e 5q cd - 
Stt 2 - e 4a 6 ab d a 6<j)d b 5(j> - e 6a V" (0)6<j) 2 - e Aa 5 ab d a 8Td b 5T 



Now, recall the Taylor expansion of a square root of two variables: 

VA + Be + Ce' = sfl- ' 



Ce' - 



(Bef 



L\~T 2\ "T V 4.1 

Thus, we can write /i p h ys as 



d 3 x Ce' 



(Bef 
AA 



Let us consider the three orders separately. 

The 0th order, , corresponds to the homogeneous Hamiltonian 



h {0) - 



G 



(108) 



(109) 



(110) 



(111) 



i.e. the Hamiltonian that generates the dynamics for the geometry in the case no perturbations are considered (FRW 
spacetime). 

As for the 1st order, notice that it involves an integral over the whole space of objects which are linear in the 
perturbations: it is not a surprise that, once we Fourier-transform it, it vanishes identically. 

The 2nd order is thus the first correction to the dynamics. After some algebra, and using the simple rules presented 
in the appendix for dealing with the Fourier transform, we obtain an object of the form 



h {2) - 

"phys - 



«) 2 



D(0)+ J2 W 
fee£-{o} 



(112) 



where D(0) contains 8q ab (0) T and 5(j)(0) T (and their conjugate momenta), while D(k) contains the k ^ modes, which 
can be expanded as Sq ab (k) T = A(k) a n b q^ n (and same for its momentum). Using the properties of all these objects, 
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and imposing the constraints and the gauge conditions, one finds that on the physical phase space r p h ys it is 



D(0) = -2Ke ia " '5 ab 5 cd 5Z ac {0y5Z bd {0y - K «\l - 5 ab S cd 5q ac (Q) T 6q bd (0y - ^5q ab (0) T 6iT ab (0y 



IS 



moyy + e ba ^'(o)w(o) 



T\2 



(113) 
(114) 



D(k)= 

m— 5,6 

l 



-2ne^( P m (kyy - ^ qm{k yp m (ky - «^ ia ( qm ( k yy - ^(<Ufc) T ) 2 



(«57f(fc) T ) 2 + e 4a ' k\8^{k) r f + e to ' V;'{0)(5j>(k) T ) 2 



(sp T {kyy - KTr T a 6T(kysp T (ky + 



(115) 



In this way, we have derived the explicit form of h p ^ ys (see (|103p ) up to the 2nd order. It is convenient to group 
the terms, according to their dependence on the fields: 



/iphys - ffhom + H k=0 + ^ H m,k + /] H k 

k=£0,m=5,6 kjtQ , 



A I 
h 



(116) 



where 

-^hom 



«(tt, 



r>2 



Hk=Q 



2 K e 4aT 6 ab 5 cd [ 6n ac (0y + * a S ae S c f5q ef (0y ) ( <y* M (0) T + ?«__s b s5 dh Sq gh (0y ) + 



12 



k(t£)V 
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<W d <5& c (0) T <%d(0) T 



t — 4a T 



12 



1 /««) 2 e 



T \2 c -4a T 



12 



4k 



(? m (fc) T ) 2 



^) 2 



l(^ T (fcr-^ r 5T(fcr 



2 + ^e AaT k 2 {ST{kyf 



+ - (e iaT k 2 + e 6 «>;'(0)) (54>(ky) 2 



(117) 



can be thought of as the various Hamiltonians generating the dynamics on the different sectors of the physical phase 
space. This concludes our exposition of the physical dynamics of the theory. 



VIII. A REMARK ABOUT THE MUKHANOV-SASAKI VARIABLES 



This is a good point to bridge our approach with the one used in the standard cosmological pertrubation theory (and, 
at the same time, showing why our is better suited in the context of quantum field theory on quantum cosmological 
spacetimc). This fact is that, if we consider only the constraints E(k), M(k), W(k), V(k) (that is the linear parts 

Cik)^, c2~\k) of the constraints of the full theory), then the gauge-fixing procedure we just presented is not necessary: 
indeed, the gauge-invariant scalar degrees of freedom are known, and are called Mukhanov-Sasaki variables [3(J. In 
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terms of our variables, they are 



W^- 2 " „(0) -2q 

Q(fc) = ST(k) + gi(fc) - gL^ g2 (fc) 



P(fc) =5p T (fc) 



((0) 6q \ 4a 

2 KTT a I KTT a 



(118) 



They are used in standard cosmological perturbation theory, because Q(k) and P(k) commute with the linearized 
constraints and have the canonical Poisson algebra: 



{Q(k), E(k')} = {Q{k),M{k')} = 0, {P(k), E(k')} = {P(k), M(k')} = 



(119) 



{Q(k),P(k')} = 8 k ,„, {Q(k),Q(k')} = {P(k),P(k')} = 



(120) 



Therefore, one does not need to fix a gauge, because the Dirac observables coordinatizing the physical phase space 
are already available: (q5(k),p 5 (k),qQ(k),p 6 (k),Q(k), P(k)), with the canonical Poisson algebra. 

However, there is a very important reason for us not to follow this route. Recall that in our approach the background 
variables a, n a , and are treated on the same footing as the perturbation variables (i.e., they are variables of the 
phase space, to be quantized). But (|1 18[) shows that Q(k) and P(k) are functions of the background variables, which 
means that whenever computing Poisson brackets involving them, one has to take into account the Poisson algebra of 
the background variables as well. Secondly, we have different constraints: namely, we also regard the contribution to 
gauge transformations coming from the second order terms of C{k) and C a (k), and we consider the constraint (7(0). 
The variables Q(k) and P(k) are not invariant with respect to all those gauge transformations. For this reason, even 
though Q(k) and P(k) are perfectly good variables for studying quantized inhomogeneous perturbations on a curved 
fixed classical background, they are not suited for the purpose of quantizing the perturbations and the background 
simultanouesly. 

Nevertheless, it is interesting to note that, with respect to these variables, the dynamics seemingly simplifies a lot. 



Indeed, let uas go back to (111 Tp . There, we grouped the terms in such a way that the Hamiltonians take the form of 
the Hamiltonian of the harmonic oscillator. More precisely: defining the new "momenta" 



5n ab (0) T := 7r a6 (0) r + 



-4a T 



12 



s ac s bd Sq cd (oy 



T — 4ct X 

p m {ky : = P m {ky + 7Ta i2 q m {ky 



(121) 



we see that (|117|) reduce to 



SP T (k) T := 8p T {k) T - -^ST{ky 



horn 



k=0 



H G - 



«) 2 



2Ke 4aT s ab s cd sn ac (oysn bd {o) 

2 K e 4QT (P m (fc) T ) 2 + 



r ( K{*lYe^ 
24 



12 



4k 



d ab s cd Sq ac (oysq bd (oy 



(q m (ky 



(122) 



Hi 



+ ^e 4aT k 2 (5T(ky) 2 



«) 2 



mm 2 i/ 4a 



+ ^ aT p + e 6a >;(o)) {6$(kyy 



From here, we see immediately how the various degrees of freedom behave: 
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• the (two polarizations of the) graviton behaves as a free relativistic particle with mass induced by the background 
geometry (via the term proportional to (7r^) 2 ) 

• the perturbations of the clock field T propagate as massless relativistic particles 

• the (perturbations of the) test field <j> propagates as massive relativistic particle, where the mass is given by the 
second derivative of its potential V$ (as it happens in flat spacetime) 

In particular, Ht is the Hamiltonian acting on the physical scalar sector, and using (|1 18[) on the reduced phase space, 
one sees that it reduces to 



6 



(123) 



In other words, if we stop at the linear order, our scalar sector (and its dynamics) coincides with the one found in 
standard cosmological perturbation theory. 

However, it is important to observe that - from the point of view of the full theory (which we are considering) - the 
transformation (|121j) is not canonical. In particular, as already pointed out, the Poisson algebra with the background 
geometry is nontrivial: 



K, m ab (oy} = ^-s ac 6 bd 6q cd (oy, {<, 6w b (oy} = * a€ * a s ac s bd sq cd (oy 

{of, p™(ky} = ^L qm ( k y, R,P m (fc) T } = ^^q m (ky 

{a\5P T {ky} = -~sf{ky, k, sP T (ky} = o 



(124) 



We thus have a dilemma: the simple form of the Hamiltonian is traded for a more complicated Poisson algebra, 
which mixes the pertubations to the clock field with the background geometry. It is important to realize this fact 
when carrying out the canonical quantization of linearized inhomogeneous modes and of the homogeneous isotropic 
background! We think that a simple kinematics is a better starting point, and thus would choose the original momenta, 
rather than the new ones (| 12 1 [) . 

Finally, notice that the test field variables S(f> T and Stt^ are real canonical variables (i.e., the only non-trivial Poisson 

brakets are {<5</>(fc) r , SiT^k'y} = 5k,k')- Therefore, for each mode of the test field <f> the Hamiltonian does have the 
canonical form of the harmonic oscillator. Due to this fact, the very result of our systematic analysis obtained in the 
current paper coincides with that of jll] obtained by using short cuts, after the restriction to 4>. 



IX. SUMMARY, CONCLUSIONS AND OUTLOOK 



In this paper, we provided a framework for quantization of linear perturbations (inhomogeneities) on a quantum 
background spacetime. We hope to have convinced the reader that, in light of canonical quantization, the classical 
Poisson algebra is the most fundamental feature to be preserved at the quantum level. This imposes a different choice 
of fundamental variables than that usually taken when quantizing perturbations on a fixed classical background 
spacetime. In particular, we showed that Mukhanov-Sasaki variables are not suited for this purpose. However, 
a natural gauge-fixing exists, which allows to use the old variables as fundamental operators, and provides a true 
dynamics in terms of the homogeneous part of the clock scalar field, T. Of course, it is expected that this formalism 
can be developed for other choices of physical time as well [3lTf33l |. 

Technically, the goal of this work was derivation of the formulae (1116j) - (|117[) for the physical Hamiltonian /i p h ys - 
Therein, the Hamiltonian is expressed by the Dirac observables 7 J (|50]) . Given value of the parameter r, the observ- 
ables parametrise the phase space of solutions to the constraints modulo the gauge transformations. The first two 
Dirac observables, a T and 7r^, are identified with the background degrees of freedom. The remaining Dirac observables 
are perturbations: 6q ab {0y , Sn ab (0y, <50(O) r , dn^Oy , 5T (k) T , 6p T (ky, q 5 (k) T ,p 5 (fc) T ,q 9 (k) T ,p 6 (fc) r , <50(fc) r , dn^ky 
where k 6 C — {0}. Their Poisson algebra is canonical (|53"| . It is defined by the proper kinematical Poisson al- 
gebra of the full theory of the gravitational field coupled to two K-G fields. The constraints are solved up to the 
first order, and the Hamiltonian itself is derived up to the second order in the perturbation variables. The physical 
Hamiltonian generates the dynamics via djj/dr = — {7J, /i p h ys }, derived up to the first order in the perturbation 
variables. The dynamics of the background degrees of freedom a T ,7r^ is an undecoupled part of the dynamics of all 
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the system parametrised by both the background variables and the variables we perturb with respect to. Going to 
higher orders in the perturbations amounts to simply adding to /i p h ys higher order terms in the perturbation variables 
and to imposing so-called linearization-stability constraints linear in the perturbations f34j - l4ll |. 

The effective difference between our results and the results of the standard approach to cosmological perturbations 
consists in the status of the Mukhanov-Sasaki variables Q(k) and P(fc): to begin with, they are not Dirac observables 
themselves in our approach, however, as any other function on the phase space they can be assigned Dirac observables 
Q(k) T and P{k) T in a r-dependent manner. They still provide the corresponding term of the physical Hamiltonian with 
the canonical form (|123[) . Nonetheless, in our framework they do not Poisson-commute with the background degrees of 
freedom (see (I124[) V This last fact has to be taken into account in the process of quantization. The consequence is that, 
whereas according to the standard approach the perturbations of the K-G field T non-vanishing in the background 
define the Hamiltonian term Ht of the same form as the Hamiltonian term Hm of the perturbations of the test K-G 
field field q> in the massless case, according to our approach the Hamiltonians take substantially different forms (see 
(|117I0 . It will make a difference between the dynamics on the quantum background of the quantum perturbations of 
the clock field T on the one hand, and the dynamics of the quantum test scalar field on the other hand. 

Our results provide a good starting point to understanding and calculating effects that quantum cosmological 
spacetime in the background has on the quantum perturbations of the metric tensor and the K-G field, specifically 
in the case of the K-G field which does not vanish in the background (in the zeroth order). 
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Appendix A: Useful formulae for the derivation of the physical Hamiltonian 

In this appendix, we collect some nontrivial formula used in the main text for the computation of the physical 
Hamiltonian. 



1. Expansion of the curvature 

The biggest trouble with (|104|) is due to the curvature terms. In this section we expand R^ 3 \ exposing S^R^ 
and §( 2 >R( 3 > in terms of the linear perturbations. The next section is dedicated to the expansion of other (less) 
troublesome quantities, namely the determinant of the spatial metric. 

The starting point is formula (I60|) 

r^ = q ab {d c v a \ d b v a \ + r a \r c c d - r/ c r 6 c d ) (Ai) 

We know that Christoffel symbols always involve derivatives of the metric, so they are at least of 1st order. Thus, the 
last two terms are of 2nd order themselves, so they are contracted by qfy, and contribute to S^R^ . On the other 

hand, the first two terms contain parts of the 1st order, so they should in general be contracted with (q?Q\ + 5( 1 * > q ab ). 

To find what the perturbation <5( 1 )g ab is, we use the definition of inverse metric, q ah qb c — &c ■ From this, it follows that 

*c° = (9(0) + 5 (1) q ab )(q^ + 6q bc ) = qfaqj® + q^Sq bc + 5^q ab q^ = 

= 5« + q$ ) 6q bc + 6^q ab q { b ° c ) (A2) 

up to 1st order. So, by contracting both sides with q^ we see that 

S (1) q ab = -QtMSqod (A3) 
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Thus, one should be not fooled in thinking that 8^q ab is simply Sq ab with its indices raised via the background metric; 
indeed, it is almost like that, except that there is a minus sign in front of it! Knowing (| A3|) . we can rewrite (|A1[) as 

R {3) = «) ~ q^)Q h ( ^qef)d c r a c b - (g$ - q^Sq^d.T a \ + q$T a d b T c c d - q$T d c T b c d (A4) 

Now, inevitably, we need to take a look at the Christoffel symbols themselves. We have 

r ab = \<f d (d a qbd + d b q a d - d d q a b) = x(gfo) + 5 {1) q cd ) (d a Sq bd + d b Sq a d - d d Sq a b) = 



2(l(o) ~ ?(o)<?(o)<%/) (daSqbd + dbSqad - ddSq ab ) 



-2a 



S cd {d a Sq bd + d b Sq ad - d d 8q ab ) - ^——6 ce S df 6q ef (d a Sq bd + d b 5q ad - d d 5q ab ) (A5) 



having used the explicit form of From this, we can compute the different objects we need: 

e -2a 

d c T a \ = — (d a d d 5q bd + d b d d 5q ad - d d d d 5q ab ) - 



(d e 6q ef ) (S d fd a 6q bd + S d fd b Sq ad - d^6q ab ) - —6q ef (S d fd a d e 6q bd + 6 d ^d b d e Sq ad - 9 e ^<5g ob ) 



p — 2a „-4a 4a 

9 b r a c c = — (F^a^gcd - — 5 ce 8 d f{d a 8qcdWqe S ) ~ —5q ef 5 ce S d fd a d b Sq cd 

e -4a e~ 4a e~ 4a 

r a d 6 r c c d = —S cd d a Sq be d e 6q cd + —6 cd d b 5q ae d e Sq cd - —6 cd d e Sq ab d e Sq cd 

-4a e _4a e~ 4a 

? d c T b c d = —5 ce 5 d fd a 5q cd d b 5q ef + —d d 5q ac d c 5q bd - —S cd d e Sq ac d e 6q bd 

(A6) 

Now, we plug these in (IA4[) . and retain only the terms up to 2nd order. The result is 

i ? (3) =5 (l) jR (3) +5 (2) jR (2) (A7) 



where 



5^R^ = e- 4a d a d d Sq ad - e- 4a S ab d e d e Sq ab (A8) 
£(2)^(3) = - 2 e- Qa 8 ab 5q bc d c d d 5q da + e- 6a S ac 8 bd 8q ab d e d B 8q cd + e~ 6a Sq ab (8 cd d a d b 8q cd )+ (A9) 

-6a 

+ e- 6a (d a 6q ab )(5 cd d b 6q cd ) ~ — (8 ab d e 6q ab )(S cd d e 8q cd )- 

p—6a o„-6a 
- —6 ab (d d Sq ac )(d c Sq bd ) - e- ea S cd (d a Sq ac )(d b Sq bd ) + -^-S ac 8 bd (d e 8q ab )(d e 8q cd ) 



2G 



2. Expansion of the determinant 

Here we present the expansion of the determinant q, also present in (|104p . One simply applies the definition (up to 
2nd order): 

3!g = e abc e de fq ad q be q cf = e abc e de f (g<°> + 8q ad )(q£> + 8q be )(q^ + 5q cf ) = 
= 3I«W + e^{8q ad $q?) + q^8q be qfj + <8q cf )+ 
+ e abc e de f(q a °j8q be Sq cf + 8q ad q^8q c} + 8q ad 8q be qf}) = 
= 3!e 6a + e 4a e abc e d ^(Sq ad S be 8 cf + 8 ad 8q be 8 cf + 8 ad 8 be 8q cf )+ 
+ e 2a e abc e def (S ad Sq be Sq cf + Sq ad S be Sq cf + 8q ad 8q be 8 c} ) = 
= 3!e 6a + 2e ia (S ad Sq ad + S be Sq be + S cf Sq cf )+ 

+ e 2a ((S be S cf - S bf S ce )8q be 8q cf + (S ad S c f - S a fS cd )Sq ad Sq cf + (8 ad 8 be - S ae S bd )Sq ad Sq be ) = 

= 3!e 6 « + 6e ia S ab Sq ab + 3e 2a (S ab Sq ab S cd Sq cd - S ab S cd Sq ac Sq bd ) (AfO) 

having used the well-known relations between e and 8. So we can write 

q = q {0} + S^q + 5<$q = e 6a + e 4a S ab Sq ab + ^l(S ab Sq ab S cd Sq cd - S ab S cd Sq ac Sq bd ) (All) 



3. Fourier mode-expansion 

Here, we explain how to get to equation (|1 12[) via real Fourier transform. Since h pi [ ys is of second order in the 
perturbations, upon plugging the expansions (|38[) in it, we will obtain something which comprises terms of the form 




(A12) 



These can be seen to reduce to 

G abcd = 8f ab (0)8g cd (0) fd 3 x + 8f ab (0)V2 £ Sg cd (k)8 kfi + 8g cd {Q)V2 J2 5 fab(k)5 kfl + 

kec + kec+ 

+ ( 5 fab{k)8gcd(k'){8 kl - k > + 8 k , k ,) - 8f ab (-k)8g cd (-k'){8 kt - k , -**,*')) ( Al3 ) 

fc,fc'6£ + 

having used the fact that 

J d 3 xe^ k - k '> x = «5 M / (A14) 

But since k, kl only take values on C+, then both 8 kt — k i and 8 k ',o vainsh when we perform the sum over k' . In other 
words, G abcd is finally reduced to 

G abcd = 8f ab (0)8g cd (0) + J2 (tfab(k)6gcd(k) + Sf ab {~k)8g cd (-k)^ (A15) 

fce£+ 

Thus, we are reduced to a sum over all k. In the case there are derivatives the result is the same: so, to expand in 
modes k, we first separate the k = mode and then simply write a sum of decoupled terms over k 6 C — {0}, each 
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of which perfectly resembles the corresponding one on coordinate space (with the difference that d a is replaced with 
k a ). s Applying these rules, we indeed recover equation (|f f 2[) in a reasonable amount of time. 
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